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The large mass limit of QCD uncovers symmetries that are not present in the QCD lagrangian. 
These symmetries have been applied to physical (finite mass) systems, such as B and D mesons. 
We explore the validity of this approximation in the 't Hooft model (two dimensional QCD in the 
large N approximation). We find that the large mass approximation is good, even at the charm 
mass, for form factors, but it breaks down for the pseudoscalar decay constant. 
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1. Introduction — Heavy mesons play a prominent role 
in our understanding of fundamental processes. From 
their weak decays we may extract fundamental param- 
eters of the standard model of electroweak interactions. 
Rare decays are sensitive to the presence of new funda- 
mental forces. And they offer the exciting possibility of 
observing for the first time violation of CP invariance in 
a decay process. 

Straightforward interpretation of measured lifetimes 
and branching fractions is marred by the difficulties that 
strong interactions present for practical calculations. 
Monte-Carlo simulations of lattice QCD may eventu- 
ally furnish accurate calculations of the matrix elements 
that are relevant to these processes. An alternative ap- 
proach is furnished by the heavy quark effective theory 
(HQET) formalism. Approximate spin and flavor sym- 
metries [1] of the ^-matrix in the one heavy hadron sector 
are made explicit. The symmetries become exact in the 
limiting case of infinitely massive heavy quarks. From 
these symmetries a number of remarkable results follow, 
such as the normalization of form factors for semilep- 
tonic B — > D and D* transitions at maximum momen- 
tum transfer, q^ ax , and a set of five relations among the 
six corresponding form factors which hold at any mo- 
mentum transfer. 

In this regard it is of paramount importance to deter- 
mine the accuracy of the large-mass approximation since 
in reality the charm and bottom quark masses are both 
only factors of a few larger than, say, the p-meson mass. 
Unfortunately, this issue involves non-perturbative ma- 
trix elements and is therefore hard to pin down. More- 
over, one has to consider each physical quantity sepa- 
rately, as the approach to the asymptotic regime of infi- 
nite masses may be faster for some than for others. For 
example, some form factors in semileptonic decay remain 
calculable at g^ ax even after 1/m corrections are in- 
cluded [2] ; and Monte-Carlo simulations of lattice QCD 
in the quenched approximation [3] indicate that 1/m cor- 
rections to the pseudoscalar decay constant of a heavy 
meson are large (of the order of 40% at the charm mass). 



In this letter we report on investigations of this issue 
using two-dimensional QCD in the 1/N expansion as a 
model of the strong interactions. It is worth empha- 
sizing from the outset that such model calculations are 
uncontrolled approximations to four-dimensional QCD, 
and as such one should refrain from using the quanti- 
tative results as estimates for physical observables. We 
believe the same comment applies to any model calcula- 
tion, like the ones afforded by potential models or QCD 
sum-rules. By gathering evidence for qualitative features 
that are common to all physically reasonable models one 
can begin to believe that such features are also present 
in QCD. 

An important limitation of this model is the lack of 
spin, and we shall have nothing to say about relations 
that follow from the spin symmetry of the HQET. On 
the other hand, because the model is a relativistic field 
theory exhibiting confinement, it is ideally suited to test 
the scale of the onset of flavor symmetries. 

In the following sections we review some salient fea- 
tures of the model, compute the pseudoscalar decay con- 
stant, compute the form factors for heavy to heavy tran- 
sitions, and close with a discussion of these results. 

2. The 't Hooft Model — We calculate the properties of 
heavy mesons in 1+1 dimensions. This model of QCD, in 
the 1/N expansion, was solved by 't Hooft [4]. It shares 
features of the four-dimensional 1/N expansion which, in 
turn, has some common ground with meson phenomenol- 
ogy. The spectrum consists of meson states, with an ap- 
proximately linear, Regge-like trajectory. Since there is 
no spin in two dimensions, these are obviously radial ex- 
citations — there is no analogue of the spin-symmetry 
relations which appear so fruitful in studies of the real 
world. Nor are there baryon or glueball states, these 
being suppressed by the 1/N expansion and the lack of 
transverse dimensions. Yet in spite of the peculiarities 
of two dimensions, it is nonetheless a non-trivial strong- 
coupling solvable theory which is ideally suited for test- 
ing our dynamical ideas. This model has been exten- 
sively studied and we refer the reader for details of the 
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calculations to [5,6] where the formalism for the matrix 
elements was derived, to [7] where the numerical meth- 
ods are discussed, and to [8], where they were recently 
applied to the study of low-energy effective theories. 

We take two flavors of quark, one heavy and one light, 
which we denote as Q and q, respectively, with bare 
masses M and to. The coupling constant g has dimen- 
sions of mass, and we work in units where g 2 N/ir = 1. 
The lowest heavy-light Qq bound state — let us call it 
the B — is a pseudoscalar of mass /x. We can solve the 
non-perturbative bound-state equation numerically for 
the wave functions and the masses [7]. 

As our interest is in the approach to asymptotia, we 
hold to fixed (taking the value m 2 = 0.3 throughout 
the paper), computing the dependence of the universal 
form factors and pseudoscalar decay constant, as 
functions of M. [10] 

3. The pseudoscalar decay constant — The pseudoscalar 
decay constant fs for the meson B is defined by 



(0\qj^ 5 Q\B(p)) = f B p» 



(3.1) 



If the states are given the usual relativistic normaliza- 
tion, 

(B(p')\B(p)) =2E5(p-p') , (3.2) 

then the large-mass limit gives the scaling behavior 
fBy/U ~ constant. This follows from the observation 
that the static properties of a 'dressed' heavy source of 
color are independent of its mass. The factor of ^/JL sim- 
ply reflects the mass-dependent normalization of states, 
cf. eq. (3.2). In four dimensions there arc logarithmic 
corrections to this relation [11]. We have omitted these 
foreseeing that they are absent in a super-renormalizablc 
theory 

In two dimensions both vector and axial currents are 
good interpolating fields for the pseudoscalar meson. 
Since -f 5 ^^ = e^ v, y v , they are both characterized by the 
same decay constant. In the 't Hooft model fs is eas- 
ily computed[5,6]. It is given by[12] fs = f$ dt(j)B(t), 
where </>s(i) is the momentum-space wave function for 
the £?-meson and t is the fraction of light-cone momen- 
tum carried by the heavy quark [4]. In fig. 1 we show 
how the limiting behavior fB\/Ji ~ constant is attained. 

Fitting the high-mass portion of the curve (M > 5) 
to a quadratic polynomial in 1/M gives a description of 
the 1/M corrections: 



f B ^ = 2 -0 



i 1.4 /1.4V 

1 -m + \m) 



(3.3) 



There are two sources of uncertainty in this calculation: 
the endpoint fit of the wavefunctions 4>b and the few val- 
ues of the mass we take for our mass fit. By varying the 
number of sampling points and successively improving 
the quality of our wavefunctions we estimate the error 



to be less than 2% for the coefficient of 1/M and less 
than 20% for that of 1/M 2 . 
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Figure 1. Plot of the pseudoscalar decay constant fByffi 
as a function of the mass of the heavy quark. The mass of 
the light quark, m, is fixed throughout: to 2 = 0.3. 

4- Form factors — Of greatest interest in probing the 
structure of the theory in the heavy-quark limit is the 
form factor for a heavy-quark current between heavy 
mesons: 

(B'(p')\Q'^Q\B(p)) = f + (q 2 Kp+ P ')»+.f-(<l 2 )(p-p')» , 

(4.1) 

where the momentum transfer is given by q = p — p' . B 
and B' have the same light quark content but different 
heavy quarks. When Q' = Q, conservation of current 
gives /_ = 0. The remaining form factor, f(q 2 ), is 
normalized, /(0)=1. 

In the limit of infinite masses, the mesons are more ap- 
propriately labelled by their velocities v — pj '/x and v' — 
p'/ri', and it is natural to describe the dependence on q 2 
through the function w = v-v' = (/i 2 + // 2 — q 2 )/2iiti! '. 
The transition amplitude for a 'dressed' heavy source of 
color with velocity v to a second one with velocity v' is 
then independent of their masses. Thus all three form 
factors in this limit are given in terms of a single univer- 
sal 'Isgur-Wise' function[l], £(W): 



f±(q 2 ) = £(W) 



Ij! ±V 



/(<7 2 )=£(W). (4.2) 



As before, we ignore logarithmic corrections [13] which 
are absent in the super-renormalizable model. 

The Isgur-Wisc function is normalized £(1) = 1 — 
by evaluating it for identical heavy mesons so that v-v' = 
1 corresponds to q 2 = 0. This in turn gives a prediction 
in the case of different heavy mesons at q 2 = q 2 lax — 

f±i&~)=(^£) ■ (4-3) 



Let us now examine the behavior of the form factors 
for finite mass heavy quarks in two dimensions, where 
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they may be computed exactly. In two and four dimen- 
sions the definitions and infinite mass relations are iden- 
tical. What are the dominant 1/M corrections and 
where do they set in? How well can the form factors 
be approximated from the quark model contributions 
alone? 

It is straightforward to evaluate the current matrix 
elements of interest. For spacelike momentum transfer 
the minus component of the current gives the following 
expression [6] in which the first term dominates our cal- 
culation: 



At 



1-t 



to 



dtt 



+ 



1 - t 

1-LU 
I - t 
l-UJ 



G(t;q 2 ) 



G(t;q 2 ) 



(4.4) 



where 

uo = p 



/<?_, the full vertex $(t) = £ At' 4>{t')/(t' - t) 2 , 
the Green function G(t;q 2 ) = f n </>QQ' (t)/(q 2 - ft), 
and G(t; q 2 ) = { At' G(t'; q 2 )/(t - t') 2 . 

The first line represents the 'quark-model' contribu- 
tion to the form factor, where the current couples di- 
rectly to the valence heavy quark. The other terms rep- 
resent the full set of remaining graphs which arise from 
gluon exchange in the current channel. These can be 
resummed explicitly and computed; we have done this 
and checked that they give but small contributions for 
the cases we studied. 

First we consider the form factor f(q 2 ) for identical 
heavy mesons, B' = B. In fig. 2 we plot the heavy-quark 
current form factors for M = 5, 14 and 450. 



Figure 2. Heavy-quark current form factors f(v-v') func- 
tions of v-v' for different masses, M , of the heavy quark. 
The solid, dashed and dashed-dotted lines correspond to 
M = 450 (effectively infinite), 14 and 5, respectively. 

To better characterize the large-mass behavior we can 
examine, at fixed w — v-v', the approach to the asymp- 



totic function as the heavy quark mass is varied away 
from M = oo. We calculated the form factors for addi- 
tional masses, M = 45 and 140. A fit to the quadratic 
polynomial in 1/M', 



f(vv') = £(W) 



1 - 



Kl(v-v') K2(v-v') 



M' 



M 



,2 



(4.5) 



gives ki = 0.10, 0.26 and 0.32, and k 2 = 0.6, 1.6 and 3.0, 
for v-v' = 1.25, 2.0 and 3.5, respectively. As was the case 
for the pseudoscalar decay constant, these quantities are 
good to about 10% accuracy. 

Next we consider transitions between different ground 
state mesons. At ^ ax the left hand side of (4.4) is pre- 
dicted in the HQET, cf. eqs. (4.1) and (4.3): 



{B'(p')\Q' 7 .Q\B(p))\ = 
In fig. 3 we show 

{B'(p')\Q'^Q\B(p))l 
d(M') = -= - 



- 1 



(4.6) 



(4.7) 



as a function of the heavy quark mass M' in the lighter 
i?'-meson. 



Figure 3. Difference d(M') from unity of the matrix el- 
ement of the flavor changing current between ground state 
mesons, at q^^, normalized to 2^7/7/7, as a function of the 
heavy quark mass M' in the lighter meson; see eq. (4.7). 



The heavier B meson has heavy quark mass held fixed 
at M = 450. This plot reveals that the finite-mass cor- 
rections are quadratic in 1/M'. This result is expected 
as a consequence of Luke's theorem[2], which states that 
there are no corrections of order 1/M' to the predicted 
normalization of form factors at maximum momentum 
transfer; corrections begin at order 1/M' 2 [14]. Again 
the result can be fit to a quadratic polynomial in 1/M': 



d(M') = - 



0.002 ±0.002 0.14 ±0.01 



M> 



M' 2 



(4.8) 



Here we have written the uncertainty as a range since 
the first term is consistent with zero. Alternatively, one 
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can fit to a power A/M' n , for which we find A = —0.17 
and n = 2.1, which clearly indicates the absence of 1/M' 
corrections. 

5. Discussion and Outlook - The size of 1/M cor- 
rections is not uniform and cannot be characterized by 
fs alone: they are far larger for the pseudoscalar decay 
constant than they are for the form factors, for which 
the asymptotic limit is approached rapidly. This is true 
both for the functional form of the form factors out to 
moderate v-v' (eq. (4.5)) as well as for the normalization 
of the form factors of flavor changing-currents (eq. (4.8)). 
A similar pattern of corrections is observed in calcu- 
lations that use the non-rclativistic potential model of 
quarks [15]. 

As we explained earlier, we do not view the model 
as a means of extracting quantitative information for 
four dimensions. Nevertheless the question of whether 
the large mass limit applies for the charm quark can be 
addressed. The transition from light to heavy quark 
dynamics was explored in detail in ref. [8]. Inspection 
of the spectrum of QQ states and the strength of the 
singularities in the form factors reveals this mass to be 
between 1 and 2 (again in units of g 2 N/n = 1) [8]. For 
this mass the 't Hooft model has corrections of order 
100% to the HQET prediction of the pseudoscalar decay 
constant, while it gives small corrections, 4-14%, to the 
normalization of form factors of flavor changing currents 
at v-v' = 1. 

It is apparent that the 't Hooft model is a valuable 
testing ground for our ideas of the large mass limit 
of QCD. One may attempt to investigate the valid- 
ity of Bjorken's sum rule [16], since one can compute 
the wavefunctions of rather high excited states of Qq 
mesons. More interestingly, one can address the ques- 
tion of whether the form factor for the semilcptonic de- 
cay of a heavy meson into a light one is wave-function 
or pole dominated. The answer remains elusive due to 
the intrinsically non-perturbative nature of dynamics in- 
volved. 

Note added — After completing this work we learned 
about similar calculations in progress by M. Burkardt 
and E. Swanson. 
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